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Abstract 

We present results of the study of the non-monotonous behavior of the Lifshitz 

^ ■ line as a function of temperature in ternary homopolymer/diblock-copolymer 

mixtures. The non-monotonous behavior of the Lifshitz line is due to the 

wave vector dependence of fluctuational corrections, which we treat in the 

framework of the renormalization group method. Our results are in agreement 

Ph ■ with the experimental findings of Schwahn et al. |l],f§. 
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I. INTRODUCTION 



> 

Q\ • The Lifshitz point appears in systems with competing tendencies for phase separation 

into bulk or spatially modulated phases. If the appropriate parameter controlling the relative 
strength of the two tendencies is varied along the critical line of phase transitions a special 
multicritical point occurs at which the character of phase separation undergoes a change from 
bulk phase separation to the phase separation into a spatially modulated phase. The Lifshitz 
point is known to exist in magnetic systems ||- ||, liquid crystals |J, polyelectrolytes 
0- m, oil/water/surfactant mixtures 0, random block-copolymers ilQ| - |TT| , mixtures of 



homopolymers and diblock copolymers [12,13]. First theoretical interest in such systems 



was generated by the work of Hornreich et al. |4j], who introduced the Lifshitz point and 
calculated the critical exponents for this class of universality. Most of the theoretical effort 
since has been concentrated on calculating the values of the exponents via application of 
various renormalization group techniques [|5],[H . 



S^ The aim of the present work is the theoretical description of the behavior of the Lif- 



shitz line with varying temperature. Recent experiments P|f2] JTo]JT^ ] on symmetric isoplethic 
A/B/A-B homopolymer/diblock-copolymer mixtures in the vicinity of the Lifshitz condi- 
tions revealed many challenging phenomena not accounted for by the existing mean-field 
theories ] |T2"| , |l"3| . In particular, one of the experimental results in clear disagreement with 



the mean-field prediction is the non-monotonous behavior of the Lifshitz line |I]J^] with tem- 
perature. We will show that the wave vector dependence of the fluctuation corrections is 
responsible for the experimentally observable shift of the Lifshitz line from its mean-field 
value. The fluctuation effects will be taken into account within the one-loop renormalization 
group method. 



We will put special emphasis on the comparison between the theoretically predicted 
behavior of the Lifshitz line and current experimental results fl],| , H| l|| . As we will demon- 
strate a major factor determining the character of this behavior is the value of the lower 
critical dimension di, a fact which has been little discussed in the literature. It is important 
to stress that the actual renormalized value of di is not known at present. Therefore, since 
the mean field value d™ = 4 is close to d — 3 — dimension of space of the considered poly- 
mer blends |I|,0,|T^,|T^] , we will theoretically analyze different types of behavior of the Lifshitz 



line resulting from different possible values of di and compare them with the experiment. 



II. PERTURBATIVE CALCULATION OF THE SHIFT OF THE LIFSHITZ LINE 

We start from the conventional expansion of the Landau free energy functional 

fT^(q)] = ^/v(-q)G f o 1 (q)V'(q) + ^/ / / ^(qiMq^faaM-qi - q 2 - q 3 ) (1) 

L Jq 4! Jqx Jq 2 J <\z 

in powers of the Fourier transform of the order parameter. In particular, for the ternary 
mixtures under consideration the natural order parameter is the deviation (from the volume 
averaged) of the concentration of a (or b, since the system is symmetric) monomers. The 
parameters of the Hamiltonian (0) for the polymer system can be obtained from coarse- 
graining of the corresponding microscopic Hamiltonian []T7| . In particular, near the Lifshitz 
line (to be defined later) the bare correlation function Gq l (q) can be written as follows: 

G 1 (q)=T + c 1 ( ( f ) )q 2 + c 2 q i (2) 

with r ~ (T—T c )/T being the reduced temperature. For the considered ternary mixtures the 
coefficient ci(0) depends on the concentration of diblock-copolymer [^2|,|13|,|l7| and changes 



the sign with the variation of 0. Within the mean-field theory the Lifshitz point is defined by 
the two conditions: r = and C\(<J)mflp) — 0. More generally, we can introduce the mean- 
field Lifshitz line as the locus of points in parameter space (<p, r) at which the quadratic term 
vanishes: c\((J)mfll(t)) = 0. The Lifshitz line (LL) is easily determined experimentally by 
considering the position of the peak of the static scattering curve [T^JITJJI|,g] . If we begin 



increasing the concentration of diblock at constant temperature the LL is determined by the 
diblock concentration at which the peak in the scattering curve first shifts off the zero wave 
vector. A noteworthy feature of the considered polymeric system is that the temperature T 
enters the Hamiltonian ([!]) only via the Flory-Huggins parameter, hence the coefficient ci(</>) 



turns out to be independent of temperature [ p?| , |T3| , |lTf , so that the mean- field position of 
the LL (J>mfll is also temperature independent and is determined solely by the ratio of the 
molecular weights of the polymers comprising the mixture. This mean-field prediction was 
not confirmed experimentally. Instead of being constant, the position of the LL was found 
to vary with temperature, more precisely, it exhibited a non-monotonous behavior, which 
shows that fluctuations should be taken into account. This is hardly surprising because at 
the LL, when ci(0) vanishes, the fluctuation corrections in fact become the only input into 
the renormalized counterpart of C\ and thus always play a role. 

Let us consider the renormalized correlation function. Note, that since our goal is to 
calculate the deviation of the LL, the renormalized quadratic term vanishes: 



G-\q) = r r + h{<t>, r)q 2 + c 2 q 4 (3) 

/ 1 (0,r) = c 1 (0) + Ac 1 (r) = O (4) 

The shift of the LL (which is temperature dependent due to fluctuation corrections) is de- 
noted by Aci(r). Within the approximation we use in this paper c 2 will not be renormalized. 
The renormalized parameter l\ can be found by considering the Dyson equation: 

G-\ q ) = G,\q)-nq), (5) 

E(q) = D x (q)+D 2 (q). (6) 

We consider in the self-energy S(g) only the one and two-loop diagrams: 



n ±lx I -1 

6 Jq r + c 2 q 
A 2 



A = -^A / -J-j (7) 



DM = T L 7 q2 [r + c 2 g 4 ][r + c 2 ql}[r + c 2 ( qi + q 2 + q)*] (8) 

For generality and ease of comparison with known results we have introduced n - - the 
number of components of the order parameter. Note, that for the polymer blends |I],[|, 
whose description is the goal of our work, due to the incompressibility condition the order 
parameter is a scalar, i.e. n — 1, as is indeed clear from the Hamiltonian ([I]). The D\ 
diagram is g-independent and is therefore of no relevance to the renormalization of C\ . The 
first correction to it is given by D 2 (q). Calculation of D 2 (q) is performed easier in the real 
space. For the experimentally relevant case d = 3 we use the r-space representation of the 
correlation function 

to rewrite the expression for D 2 (q) in terms of G(r) as follows: 

D 2 (q) = ^Je^G 3 (r)d 3 r (10) 

where £ = (c 2 /t) 1 ' 4: is the mean-field correlation length. In fact we need only the quadratic 
term of the diagram in powers of q, which is readily calculated: 

Df\q) = -b^ (11) 

c 2 

where b ~ 0.109 10~ 4 is a constant. As it is clear from the Dyson equation (|5|) this result 
gives in fact the shift of the LL: 

\2cS 

Aci (r) = 5-^ (12) 

C 2 



Note, that we have obtained the expression ([12]) within the perturbation theory and therefore 
it is not valid in the regime of strong fluctuations. However, the scaling behavior of Aci(r) in 



the regime of strong fluctuations can be obtained from (|i"2|) by replacing A with the effective 
coupling constant X r , and understanding under £ the true correlation length, thus 

AC1(T) . B *)|!M. (13) 

C 2 

This formula immediately allows some conclusions about the qualitative behavior of the LL. 
First of all, we observe that the correction is positive, which means that on the LL Ci(0) < 
0, i.e. fluctuations shift the LL into the g* 7^ region of the mean-field theory. For the 
homopolymer/diblock copolymer blend that means that LL shifts to greater concentration 
of diblock, which is in agrement with experiments |15| , [16| ,[1|,|^| . Next consider the dependence 
on temperature. We have two regimes here: perturbative (small correlation lengths, at high 
temperatures) and scaling (low temperatures, large £). In the perturbative regime where the 
input of fluctuations is small, A remains practically non-renormalized, so that with lowering 
temperature Aci(r) should increase simply due to the increase of £. In the scaling regime the 
main effect (as will be shown below) comes from the renormalization of the coupling constant 
A r . In fact in this regime we can obtain the scaling dependence of the correction from the 
considerations of dimensionality: demanding that the Aci(r)g 2 term of the Hamiltonian ([3]) 
have the same dimensionality in £ as the C2(f term of the correlation function. If c<i is not 
renormalized (as is our case), then 

A Cl (r)~r 2 (14) 

As we can see, in the scaling regime the correction decreases with increasing £ (i.e. decreasing 
temperature). Combined with the conclusion made above about the increase of Aci(r) in the 
perturbative regime we come to conclusion that Ac\ (r) exhibits a non- monotonous behavior 
as a function of r. This behavior is a manifestation of the crossover between the regimes of 
small and strong fluctuations. 



III. RENORMALIZATION GROUP STUDY OF THE LIFSHITZ LINE 

To describe Aci(r) quantitatively we have to obtain expressions for A r (r) and £(r) in both 
the perturbative and scaling regimes. For this purpose we shall employ a renormalization 
group technique to the first order in e (one-loop RG). Note, that within this method the 
parameter c 2 does not renormalize. The renormalization of temperature is described by the 
one-loop diagram (|7|): 

D 1 = Dt + D[ = -?±1\ f^- 2 + ^r\ J - \ (15) 

6 iq C 2 (q^ 6 JqC2(q) (1" + C2(q)) 

The above expression is conventionally split into two parts responsible for additive and 
multiplicative renormalization of temperature. For dimensions d > 4 a cutoff at the upper 
limit in integration over q in the first term is assumed. These two terms give the critical 
dimensions of the Lifshitz class of universality. The lower critical dimension d\ is defined as 
the dimension when the first term logarithmically diverges at small q. The upper critical 
dimension d u is the dimension at which the second term logarithmically diverges at small q 
for zero temperature. A cutoff at the lower limit of integration over q is implied in Equation 



(|l5l). For the isotropic Lifshitz class of universality we obtain: dj = 4 and d™f = 8. 
The real experimental system corresponds to d = 3 so that we come to conclusion that we 
are situated below d™ . This means that the D^ term diverges at small q, i.e. for large 
correlation lengths and thus no phase transition of the second order is possible at a finite 
temperature. However, this value for the lower critical dimension is only the mean-field one. 
Fluctuations renormalize the value of the lower critical dimension. The calculation of the 
renormalized lower critical dimension is a formidable task, so that the true renormalized 
value of di for the experimental system is actually unknown. Therefore, we will consider 
below several possibilities. 

Going back to renormalization of r by substituting ([15]) into the Dyson equation we 
obtain: 

T r = T a Z 2 (A min ) (16) 

r a = T-D a 1 Z 2 (A miQ )- 1 , (17) 

where for purposes of clarity we introduced the temperature with additive term r a as well 



as renormalized temperature r r . The quantity A m j n in Equation (|16|-|T7|) is the lower cutoff 



imposed in Equation (15) in integration over the momentum q. The RG treatment is based 



on the following perturbative expression: 

1 




q c 2 (q 2 ) 2 (r + c 2 (q 2 ) 2 ; 



+ ... (18) 



Likewise considering the fluctuation correction to A for its renormalized counterpart A r we 
obtain: 

A r = A (l - ^A I V^-m + -) ( 19 ) 

V 6 V q (r + c 2 (q 2 ) 2 ) 2 J 

These two equations are the starting point to derive the differential equations of the renor- 
malization group for r r and A r . To do this we introduce a running cutoff A at the lower 
limit of the integrals in ([18| ) - ([T9|) , differentiate both parts of these equations with respect to 
A and replace in the rhs the bare quantities r and A through the effective ones. Thus we 
obtain: 

, d In r r n + 2 , . 

dA 6 y v ; 

A_3 = -eg + g 2 21 

dA y y 6 y K ' 

where the effective dimensionless coupling constant is defined as follows g = A r A~ £ with e = 
8 — d, and A = XKa/cl, fQ = Sd/(2n) d , Sd being the surface of a unit <i-dimensional sphere. 
As can be seen from (|2T| ) the fixed point of the effective coupling constant, g = A r A~f n , is 

9* = 4^- 

v n+8 

Notice that the second equation is independent of r, therefore we solve it first and then 
substitute the result g(A) into the first one. Thus we obtain the solution of the differential 
equations ©-©: 



'a 




j = Z 2 (A„i„) = (l + ^^KL ] (22) 



^( 1 + ^ A3 „)-* (23) 

Equation ( p2| ) for r r allows us to obtain the critical exponent of the correlation length: 
v = | f 1 + ^| £ ) , which is a well-known result @|J. Note, that in this system r r is expressed 
via the temperature with additive shift (RG generalization of Equation (|ITD): 

6 Jq Z 2 {q)c 2 q^ 

Taking into account the renormalization of the coupling constant A r and the propagator 
in the expression of the shift of the critical temperature can be found by considering the 
higher-order corrections to the self-energy £(<?). It is clear that the infrared behavior of 
these corrections is controlled by the momentum q, which is the argument of the self-energy 



S(g). This demands to write X r (q) and Z 2 (q) under the integral in (|24j) as functions of the 
external momentum q. The relation ( |2~3| ) makes the one-loop RG scheme for renormaliza- 
tion of r r and A r complete. Now, using the relation between A min and £ (see below) we 



can obtain from Equation fl22[) £(r), which substituted into (|23|) will give A r (r). The two 



dependences substituted in turn into the formula for Aci ([Uf) will give our final result - 
the deviation of the LL from the mean-field value as a function of temperature. In order to 
find relation between the cutoff wave vector A min and the correlation length £, one should 
find the perturbative limit of the RG formulae (|22|) or ( P3|) and demand it to be equal to 
the corresponding diagrams ( p~8| ) or ([19]). Thus it is straightforward to obtain: A min = £ _1 . 
Using this relation and introducing reduced variables we can rewrite Equation fl22l) -(f24|) as 
follows: 

5 = z 2 (0 = (i + xeY m , f a = f + a r ^rjq'-'dq (25) 



Tn 



r 1 Z 2 (q)q* 



y = (l + Ar)" (26) 

where the reduced variables are: A = lt ^--, f = — , and the constant a = ^rfe. However 
the relation A min = £ _1 is asymptotically correct only in the vicinity of the upper critical 
dimension d u = 8 and we do not expect it to hold for the considered experimental situation 
d = 3. Therefore, we only know that A m j n ~ £ -1 with the prefactor being unknown. In this 
situation the constants a and a new constant / (defined via A = /A/c|) become essentially 
fit parameters of the theory. The expression for Aci (r) in reduced variables reads: 

Ac l( r ) T2/~x t 8 



c 2 bt 



K(r)e(r) (27) 



Before solving the system let us make some further comments on the Equation (p2|)-(p4[), 



in particular discuss the issue of the lower critical dimension. To that end we should consider 
the scaling of the term responsible for additive renormalization: 

6 



^ - « /" 7 (\u \ * qd ~ ldq ~ f 17 " < 28 ) 

If d > di then this integral converges on small g and this addition to temperature can 
be neglected in the RG calculations as it is does not depend on the correlation length. 
Note, however that experimentally this correction is still relevant. Since the values of the 
fluctuation shift of temperature are different in the Ising and Lifshitz classes of universality 
the temperature of the LP is shifted to lower values then the transition temperature of 
the Ising class of universality. The situation of the lower critical dimension corresponds to 
logarithmic divergence of At, i.e. 1/v = 0. If d < di then the correction diverges at £ — > oo, 
which precludes the phase transition at finite temperatures, the transition temperature goes 
to zero. According to Equation fl28|) we can consider the exponent v in this case to be 
formally negative. As we have mentioned above the renormalized value of d\ is not presently 
known, so we will consider the two possibilities: d > d\ and d < d%. Therefore, it does 
not make sense to use the exponent v obtained to order e, since we do not expect it to 
be correct for d = 3. On the contrary, using the scaling relations (supposed to be correct 
even for negative v) we will express the exponents in the equation for £(f) fl25|) via the 
exponent v (generally the exponent rj is also necessary, however it is zero in the one-loop 
approximation). Thus, we arrive at the equation 

(r + aXlZ (1 + A^)- 1+<4 - 1M/ * t^y =?(! + Ar)- ,4 - 1/ "* / * (29) 

Its solution £(r) substituted into ( PS| ) gives A r (r) and thus we can obtain the shift of the 
LL Aci(f) according to (p7|) . 

The results of numerical evaluation of Aci(r) according to Equation (p7| - p9| ) are plotted 
in Figure 1. We have considered two cases: d > di, v = 1 — solid line (A = 1, o = 2); d < d\, 
v = — 1 — two dashed curves (A = 1 for both curves; a = 0.5, a = 2). If d > di then for all 
values of parameters a and A the LL has the qualitative form as the curve plotted in Figure 
1: with decreasing temperature Aci(f) initially increases (due to increase of £ as explained 
above), then for large fluctuations decreases (due to decrease of the renormalized coupling 
constant A r ) and the LL ends at a Lifshitz point. Note that the LP occurs at the negative 
temperatures f due to the finite shift of transition temperature discussed above. For d < di 
the situation is more complex. At large values of a (a = 2 curve of Figure 1) the LL goes 
tor^ — oo, so that no Lifshitz point exists (we may say it is shifted to f = — oo). As can 
be seen from Figure 1 in other respects the LL of this case is qualitatively the same as in 
the case d > d\. At small values of a (a = 0.5 curve of Figure 1) two solutions of Equation 
( p9|) exist at large temperatures (the one corresponding to greater £ is plotted with the dot 
curve in Figure 1) and no solution at small temperatures. Thus in this case the LL ends 
in a point at which the correlation length is finite, which indicates that this is a point of 
the phase transition of the first order. It is not clear however, whether this case represents 
physical features or is an artifact of approximations of the theory. 

Comparing obtained results with the experimental data of Schwahn et al. |l]|| we observe 
that the experimental LL has the same qualitative shape. We would like to stress that the 
present study predicts that the Lifshitz line approaches its mean-field counterpart for high 
and deep temperatures, the prediction which is also in agreement with experimental behavior 



of the Lifshitz line found in |],|2|]. However, on the basis of the current data of Ref. |l],|2j it 
is not clear which of three cases concerning the value of the lower critical dimension of the 
system which we have discussed above, actually takes place. 

Summarizing, we have calculated the deviation of the LL from the mean-field behavior by 
taking into account the first fluctuational correction to the C\ coefficient. The renormalized 
properties are calculated within the one-loop renormalization group. Our calculations are 
carried out in the immediate vicinity of the Lifshitz line, so that no crossover to either Ising 
or Brazovsky universality class is considered. The obtained behavior of the LL qualitatively 
agrees with that observed in the scattering experiments on ternary homopolymer/diblock- 
copolymer systems. The approach we have presented in the present article is based on 
the Landau- Ginzburg expansion of the free energy with phenomenological parameters A, 
Ci, C2- Of course, this fact restricts a complete quantitative comparison with experiment. 
However, we point out that despite this the agreement of the behavior of the computed 
Lifshitz line at high and low temperatures with experimental one is a strong support of the 
validity of our approach. The major reason for the uncertainty in the complete quantitative 
comparison with the experiment is due to unknown values of the critical exponents of the 
Lifshitz universality class. Depending on the value of the lower critical dimension in the 
renormalized theory the LP exists or it is shifted to infinitely small temperatures. More 
experiments are needed to determine which of the possibilities actually takes place, which 
would amount to experimental determination of the lower critical dimension of the system. 
We hope that this work will stimulate such investigations. 

The authors acknowledge stimulating discussions with E. Straube, I. Erukhimovich, and 
D. Schwahn. A. K. acknowledges support of the Graduirtenkolleg "Polymerwissenschaften" . 
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FIGURE CAPTION 

Figure 1. The shift of the Lifshitz line Aci as a function of temperature f. Solid curve: 
d > di, dash curves: d < d L . For all curves A = 1. 
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